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We develop a theory to describe the damping of elementary excitations of a Bose-condensed gas in 
the hydrodynamic regime for the thermal cloud. We discuss second sound in a spatially homogeneous 
gas and the lowest excitations of a trapped condensate. 



The discovery of Bose-Einstein condensation (BEG) in 
trapped alkali atoms vapors has stimulated a grow- 
ing activity in the field of ultra-cold gases Q. Recent 
developments allow one to reach sufficiently high densi- 
ties and temperatures of a Bose-condensed gas, where the 
thermal cloud leaves the collisionless domain and enters 
the collision dominated (hydrodynamic) regime Q . This 
is the regime that has been studied in superfluid liquid 
helium. However, due to the high density and the com- 
plexity of interparticle interaction in helium, only few 
results can be fully described on the basis of an ab initio 
theory. In this respect, the study of hydrodynamic ex- 
citations of a Bose-condensed gas is important, since it 
bridges together phenomenological models of condensed 
matter physics and microscopic theories. 

The damping of elementary excitations of a Bose- 
condensed gas in the hydrodynamic regime is qualita- 
tively different from that in the collisionless regime. In 
the latter case the excitation frequency u; greatly exceeds 
a characteristic relaxation rate in the thermal cloud , 
and this cloud behaves as an ensemble of Bogolyubov 
quasiparticles. Theories of damping phenomena for a 
Bose-condensed gas in the collisionless regime are re- 
viewed in Q and for the earlier work in Q . 

In the hydrodynamic regime, where ivth <^ 1, the ther- 
mal cloud rapidly responds to the motion of the conden- 
sate. Accordingly, quasiparticles of the cloud loose their 
individual character, and the spectrum changes. Then 
the dynamics of a Bose-condensed gas is described by 
the Landau-Khalatnikov (LKh) two-fluid hydrodynam- 
ics characterized by two independent branches of el- 
ementary excitations: first and second sound. Recently 
a microscopic derivation of the two-fluid hydrodynamics 
has been generalized for the case of trapped gases (see ^] 
and refs. therein). For a homogeneous Bose-condensed 
gas a microscopic theory equivalent to the LKh hydrody- 
namics was developed by Popov . Thus far this theory 
has only been used for calculating transport coefficients 
at very low temperatures . 

In this Letter we develop a theory to describe the 
damping of low-energy excitations of a Bose-condensed 
gas in the high-temperature hydrodynamic regime. We 
discuss second sound in a spatially homogeneous gas and 
the lowest excitations of a trapped condensate at temper- 
atures T ^ fi (the chemical potential fi = griQ, with hq 



being the condensate density, g = ■iirh'^a/m, m the atom 
mass, and a the scattering length) and assume that the 
density of non-condensed atoms n' <^ hq. In this case, 
second sound representing out-of-phase oscillations of the 
condensate and thermal component, is mostly seen as 
the condensate oscillations, since the amplitude of oscil- 
lations of the thermal component is much smaller (see 
|lO|). The friction between the condensate and the nor- 
mal component, originating from collisions between Bose- 
condensed and thermal particles, leads to energy transfer 
from the condensate to the thermal cloud and, hence, to 
damping of the condensate oscillations. 

We will start with calculating the damping of second 
sound in a homogeneous gas. At T ^ /i and n' <^ hq 
the velocity of second sound is eg = \/nog (hereafter 
h = m = 1) and it is much smaller than the velocity of 
first sound u ^ \/T ^ representing oscillations of the 
thermal component. In this sense, the two branches of 
the excitation spectrum are adiabatically separated. For 
these conditions, the deviation of the condensate den- 
sity (5rio(t,r) from its equilibrium value uq satisfies the 
hydrodynamic equations [ pd| 

'-^^+noV. = 0, m'^+gVSnoit,r)^0, (1) 

where v is the superfluid velocity. The solutions of Eq. (|^) 
are sound waves with momentum k and frequency uj^ 
satisfying the dispersion relation Uk = csk. 

As we will see below, the damping of the condensate os- 
cillations is mostly provided by thermal excitations with 
energies e T ^ fj,. In a dilute gas they are always in 
the collisionless regime and, hence, the excitation energy 
is ep = + gno{t, r), where p is the momentum of the 
excitation, and no(i,r) = uq + Sn(){r,t) is the slowly os- 
cillating condensate density. We will assume that these 
thermal excitations are quasiclassical and describe them 
by using a classical distribution function f(t,r,p). 

The behavior of the thermal excitations is governed by 
the Hamiltonian Hq + -f^int' where 
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Oscillations of the condensate density can be written as ((1 + /i)/2/3 — /i(l + /2)(1 + fs)) + S{fp^ — ep^ + ep j 



Sno{t, r) = i 



k 



1/2 



(cke'^"""'""* - c-c), (4) 
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-P3((l + /l)(l + /3)/2~(l + /2)/l/3)], 



with Ck(c^) being a classical limit of the cre- 
ation(annihilation) operators of second sound phonons. 
It is the interaction term -ffjnt® ^^^^ leads to energy 
transfer from the condensate oscillations (phonons of sec- 
ond sound) to the thermal cloud. 

We will perform our analysis within a linear response 
theory. In the presence of phonons, the equilibrium 
distribution function of thermal excitations /o(p) = 
(exp(ep/T) — 1)~^ acquires a small extra term 5f. The 
Fourier transforms of 5f{t, r, p) and (5rio(t, r), 



where fi = f{t,ri,Pj), and the factor 1/2 in front of the 
first term in the rhs appears due to the double count- 
ing of final states. In the spatially homogeneous case 
the condition n' <C no requires temperatures well be- 
low the BEG transition temperature Tc, whereas in a 
trapped gas it is satisfied even at T close to Tc (see, 
e.g., Q]). Since in the linear regime the deviation 5f of 
the distribution function from the equilibrium value /o is 
small, we linearize the collisional integral /p. Introduc- 
ing the function hp{t,r) = 6f{t,r,p){dfo/dep)~^ and 
setting /p[/o + hdfo/dep] = ~{dfo/dep)ip[h], we obtain 



( 7not:^) ) = / '^^^^^^P - ( Yi!(''? ) ' /p, [h] = ng^noil - exp (-epjT)) J] /o(.pJ/o(.p 
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are related to each other by means of a generalized sus- 
ceptibility ap(ti;, k) fl^ : 

<5/p(w,k) = ap{Lj,k){dfo/dep)gSno{uj,k). (5) 

The imaginary part of ap(w, k) is responsible for dissi- 
pation and is related to the damping rate Tk of the (first 
sound) phonon with momentum k and frequency LUk : 



Tk = Wfc.glm((ap(a;fe,k))). 



(6) 



Eq.(||) follows from the fact that the energy loss rate 
of the condensate oscillations is given by X^k '^'^I'^kl^^'^ 
and should be equal to the rate of the energy in- 
crease in the thermal cloud. The latter reduces to 
H = Y.k^k\gSno{uj,k)\Hm{{ap{uj,k))), where ((...)) = 
Epdfo/dep{...) (see 0). 

The response function ap can be found from the solu- 
tion of the kinetic equation for the thermal excitations. 
As the force acting on a thermal particle is — (^VriQ, the 
kinetic equation takes the form 



df df 



df 



dr 5Vno^=/p[/], 



(7) 



where Ip[f] is the collisional integral. There are two 
types of binary collisions contributing to Ip [f] : collisions 
of non-condensed with Bose-condensed atoms, transfer- 
ring the latter to the thermal cloud (and vice versa) , and 
collisions between non-condensed atoms, leaving both of 
them in the thermal cloud. The former process is char- 
acterized by the relaxation rate ^ noavx, where 
a = Stto^ is the elastic cross section and vt the ther- 
mal velocity, and the latter process by the rate n'avx- 
Assuming the inequality n' <^ no, we confine ourselves 
only to collisions involving Bose-condensed atoms. The 
corresponding collisional integral is 



Ip^ [/] = 2^g2^o E[(l/2)'5(ep, - ep, - ^pJSp,,p,+Ps 
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X 



(exp (ep^ / T)6{ep^~ip~ ep^ )Sp^ ,p^+p^ {hp- hp- hpj~ 



2 exp(ep^ /T)6{ep- ep^+ (pJSp^,p-p^ {hp- hp^+ hp J). 

(8) 

The collisional integral has a number of important 
properties 0. ColHsions conserve the total momentum 
and energy of particles, and thus /[p] = and I[ep] = 0. 
In addition, the collisional integral is orthogonal to the 
integrals of motion: For an arbitrary function x we have 



(((f KpM» = ([J). 



(9) 



We now use Eq.(|^) and rewrite the kinetic equation 
(0) in the Fourier representation: 



L[ap] = {u - kp)ap + kp = -iip[ap]. 



(10) 



Eqs. (|To|) and (|) allow us to find the response of the ther- 
mal cloud to the condensate oscillations and to calculate 
the damping rate for the modes of second sound. In the 
hydrodynamic regime {lotji <C 1), Eq.([l0|) is dominated 
by the collisional integral and can be solved by using the 
Chapman-Enskog method (we are following the Popov 
ansatz for a Bose-condensed gas Q). We will seek the 
solution of Eq. ( |lo| ) in the form of expansion in powers of 
a small hydrodynamic parameter luth: 



- JO) 



v(l) 



(2) 



(11) 



o^Tfl ^ 1. The collisional integral 



where ap jctp 
can be estimated as I\a\ ^ ol/th and, hence, we have 
/[ofp'*'^''] ^ LOOp . Therefore, Eq.([l^) is equivalent to an 
infinite set of coupled equations: 



2 



(0)i 

P J 



(12) 
(13) 



In the so called acoustic approximation, one neglects 
all terms with i > in Eg. (pi]). The zero order solution 
afi follows from Eq.(|l|) and is a linear combination of 
the integrals of motion: 



(0) 



6kp. 



(14) 



The functions a and h should be found from the substitu- 
tion of Eq. ([l^) to the orthogonality conditions (^ . This 
gives a system of two linear inhomogeneous equations: 



where (\)\ and (/>2 are solutions of the equations 
/p[(p.P.-f/)'/'i] = te-^/), 



4>(P'/'2 



(18) 
(19) 



The functions (5a and bh are obtained from the require- 
ment that both sides of Eq.(13) (with i — are orthogo- 
nal to the integrals of motion. The corresponding system 
of equations is similar to Eqs.(p^,(|T^) and gives 



(5a(tjfe,k) ~ 



6u^ V ((P'» 



-4- 



2 ((pVi)) 

5 m) 



(20) 



(15) 
(16) 



Eqs. ( p^ ) and (^6|) describe the propagation of fluctu- 
ations in the hydrodynamic thermal cloud. Without the 
rhs, the solutions of these equations are free oscillations 
of a gas of excitations (first sound). They are charac- 
terized by the dispersion relation = fc^u^, where the 
velocity u of first sound is given by 



In the temperature domain T 3> /i we have = 0.86T. 

For finding the damping rates of first sound from 
Eq.(^ we are interested in frequencies lo ^ ujk ~ csk 
which at T 3> /i are much smaller than = uk. For 
these uj the solutions of inhomogeneous Eqs.(p^,(p^) are 
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(17) 



and from Eq.(p4|) one gets a real response function. Thus, 
in the acoustic approach, Eq.(^) gives zero damping rate. 

The damping of second sound appears in the so called 
viscous approximation accounting for the first order cor- 
rection ap"*. The latter is a solution of Eq.(|l^), with 



and aj?^ from Eqs. (|T|), (0) 



kp(ijj6 — aep 



1) - ^ bkikj{piPj 



p2)=^/p[aW]. 



As the collisional integral vanishes for ap ■' cx ep and cx 
kp, the solution of this equation can be written as 

a'p = Sap + Saep + ^fo(kp). 
The quantity Sap can be represented in the form 

Sap = z(kp)(/)2 -t- iuipi 
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Both Sa, Sb and 0i, (t>2 turn out to be independent of the 
direction of the vector p. Therefore, one easily finds that 
only the term Saep contributes to the quantity {{ap^)). 
Then, using Eq.(|^), for the damping rate of second sound 
we obtain Tk = U!kgini{{Sa(ujk,'k.)ep)) . To calculate the 
averages containing 0i and 02 in Eq.(pO|), we invert nu- 
merically the collisional integral in Eqs. ( p^ ) and (p^). 
Then, at temperatures T 3> /i = nog we find 



4.8—^. 

nog 



(21) 



The main contribution to (5a(wfc,k) in Eq.(|2(]|) and to 
the damping rate Tk ( pl| ) is provided by thermal exci- 
tations with energies ~ T. As the characteristic relax- 
ation time for these single-particle excitations is t/} ~ 
{nocrvT)^^ , the criterion of the hydrodynamic regime for 
the damping of phonons with frequency tOk takes the form 



(22) 



The hydrodynamic criterion (|2^) is not exactly oppo- 
site to the condition of the coUisionless regime, cuk ^ 
T{nQa^y/^ , found in [n5|. In the latter case 



the collisional integral in Eq.(hOf) is small and ap 
-kp/(w-kp). Then, Eq.(|) gives the damping rate 



r 



T 



nog 



(23) 



which, except for the numerical coefficient, coincides with 
the rate of Landau damping in the coUisionless theory 
p5[ . The main contribution to F^,^'' ( p3| ) comes from the 
interaction of the oscillating condensate with thermal ex- 
citations which have energies ^ n^g [|16| and are charac- 



terized by the relaxation rate 
exceeds the relaxation rate 



For T » nog it 
of atoms with energies 
(nQglTyl"^ <^ 1. Hence, there is a range 



of frequencies, where ^ W/t ^ t^^^^ and neither the 
hydrodynamic nor the coUisionless approach is valid. 

The comparison of Eq. (|2^) with Eq. (^) shows that the 
ratio of the hydrodynamic to coUisionless damping rate 
is Ffc/F^,^-* ~ (jJkTR, i.e. in the hydrodynamic regime the 
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damping is much slower than one would estimate with a 
collisionless approach. 

We now calculate the hydrodynamic damping of the 
lowest excitations of a harmonically trapped Thomas- 
Fermi condensate. In this case, for the density fluctu- 
ations we have (see El) 



Sno{t, r) = i 



1/2 



W,{r){c,e' 



c.c.}. (24) 



For a low-energy (ui, <C /i) excitation with a set of quan- 
tum numbers i', the function W^, is a polynomial of the 
coordinate variables ri. The condensate density has the 
well-known parabolic shape no(r) = nom(l — 
with the maximum value nom = 7^o(0) = fj./g and the 
size k = in the i-th direction. 

The wavelength of the lowest excitations is of order 
the spatial size of the condensate, which in the hydro- 
dynamic regime greatly exceeds the mean free path of 
thermal excitations. Therefore, we will use a local den- 
sity approach assuming that in the vicinity of a point r 
the generalized susceptibility is the same as in a homo- 
geneous gas with the condensate density no(r). This is 
justified by a local spatial character of the collisional in- 
tegral. As Eq.( |24|) only contains an extra factor Wy(r) 
compared to EqT(Q), the damping rate is given by r,y = 
Lv^g J dr|W^jy(r)pIm((ap(a;j/, k, no(r)))) and reduces to 



(25) 



r, = / dr|W^,(r)|^r,;,(r). 



The quantity T^hi^) is the damping rate of first sound 
phonons with frequency in a homogeneous gas with 
the condensate density no(r). This quantity is given by 
Eq.(pTl) with no — no{r), and the integral in Eq.(|25|) 
is logarithmically divergent at the condensate boundary. 
The divergency originates from our assumption that uq is 
much larger than the density of non-condensed atoms n' , 
which is violated near the Thomas-Fermi border. Hence, 
we have to cut-off the integration at distances where 



Then, expressing Wi, in terms of the ra- 



dial variable y — (^irf /IfY^"^ and angle variables 0,(p, 
with logarithmic accuracy we obtain 



2.47,. log! — r ), 



(26) 



where 7^ = / dVt\Wu{y = 1, 0, (^)|2. 

Our approach describes the damping of oscillations of 
a trapped condensate, assuming that the latter was set 
into motion while not directly affecting the thermal cloud 
(or vice versa). Using the polynomials from [ p^ , we 
calculated the coefficients 7^ for several particular cases. 
In spherical traps we find 7^ — (4n -1-2^-1-3), where n 
is the radial quantum number and I the orbital angular 
momentum. For the so called surface excitations (n = 0) 
this result remains valid in cylindrical traps, and Eq.(|2f 



indicates a significant increase of the damping rate with 
increasing /. However, it is less dramatic than observed 
at MIT in the collisionless regime. For the dipole 
axial oscillations of the condensate we obtain 7^ = 5. 
Then, for the conditions of the MIT experiment [|j {lOz = 
18Hz, T « l^K, ^l « 200nK), Eq.(||) gives F^ « 10 s^^. 
This agrees with the experiment which was entering the 
hydrodynamic regime for the axial motion. 

We believe that our approach can be further developed 
to describe damping phenomena in various intermediate 
regimes which one can be meet in trapped gases. 
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